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We investigate the role of pressure in a class of generalised Skyrme models. We introduce 
pressure as the trace of the spatial part of the energy-momentum tensor and show that it 
obeys the usual thermodynamical relation. Then, we compute analytically the mean-field 
equation of state in the high and medium pressure regimes by applying topological bounds on 
compact domains. The equation of state is further investigated numerically for the charge one 
skyrmions. We identify which term in a generalised Skyrme model is responsible for which 
part in the equation of state. Further, we compare our findings with the corresponding 
results in the Walecka model. 

PACS numbers: 


I. INTRODUCTION 

The Skyrme model [1] is a low energy effective model of QCD where baryons (atomic nuclei) are 
realised as emergent objects in a “mesonic fluid”. In modern language, they are solitons in a held 
theory based entirely on some mesonic held U, in the simplest version, U G SU{2). The pertinent 
topological charge carried by a soliton (skyrmion) may be identihed with the baryon charge. 

The most general Lagrange density respecting Poincare invariance and leading to a standard Hamil¬ 
tonian reads 


-^0246 — -^0 + -^2 + -^4 + -^6 ) (I-l) 

where 

Ta =-AsTr £4 = A4Tr ([L^, Ce =-{247r^fXGl3^I3^^. (1.2) 

Here, the left invariant current is given by 

= U^d^u (1.3) 

and Co is a non-derivative part, i.e., a potential. In its simplest version £ 024 , where the sextic part 
is neglected, the model was shown to describe baryons and light nuclei with a very good accuracy 
m-m- However, this particular proposal leads to results which are difficult to reconcile with two 
important qualitative properties of atomic nuclei and nuclear matter. It gives rather large binding 
energies and describes crystal-like nuclear matter [8j, whereas nuclear matter is only weakly bound 
in atomic nuclei and behaves more like a huid. These unwanted properties of the usual Skyrme 
model can be resolved if one adds a larger amount of the sextic term. This is related to the 
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observation that there exists a submodel within the family of Skyrme type models which has the 
BPS property (classical zero binding energies |^) and describes a fluid (has the energy-momentum 
tensor of a perfect fluid [lOj and is, in fact, a field theoretical realisation of the Eulerian fluid HU)- 
This submodel is called the BPS Skyrme model [9j and reads 

k^BPS = C&+ £-0 (1-4) 

where Cq is a further potential. Then, an improved proposal which might also be valid for the high 
baryon charge regime, is given by the near-BPS Skyrme model | 12 j 

C = e (To -|- T 2 + T 4 ) -|- P-BPSt (1-5) 

where the usual Skyrme model only provides rather small corrections to masses (binding energies), 
as e is assumed to be a small parameter. Let us remark that the sextic term Cq can be effectively 
introduced by a coupling with the w meson |13j-|18]. Here we consider the general Skyrme theory 
with the SU ( 2 ) fields only, where some effects induced by further mesons emerge due to a particular 
form of the action rather than by the inclusion of some new explicit degrees of freedom. 

Although the BPS limit is physically very well motivated and provides rather accurate results 
for binding energies for larger nuclei (after a careful treatment of semiclassical quantisation of 
(iso)-rotational modes, inclusion of the Coulomb interaction and iso-spin breaking [ 12 ] (see also 
Ref. [19])), it is still necessary to investigate the full near-BPS version. The inclusion of the non- 
BPS part leads to the appearance of physical pions and determines the proper geometry (shape) of 
skyrmions [ 20 ], which is of high importance as it results in particular patterns for the iso-rotational 
excitations. Hopefully, one can get a model which unifies the convincing results for baryons and 
light nuclei with the crucial properties of the BPS model. This is, in principle, a difficult task since 
the solvability property of the BPS model is lost once the perturbative part (usual Skyrme action) 
is added. Moreover, such a perturbation is non-analytical as the non-BPS part is the dominating 
part close to the vacuum. Hence, no small e expansion exists at the level of the field equations. 
Some recent findings for solitons with first few baryon charges, however in a not so “near-BPS” 
regime, can be found in Ref. m (see also [ 22 ]). 

The aim of the present work is to investigate the issue of the equation of state (EoS) in the 
general Skyrme model with a special focus on the near-BPS model. The BPS submodel allows for a 
rather complete analysis of its thermodynamics at zero temperature, as a consequence of the BPS 
property and of its perfect fluid form which, however, no longer hold for the general (non-BPS) 
model. The first difficulty, therefore, is to properly define the pressure, since in the absence of a 
perfect fluid description no obvious definition exists. This is analysed in section 2. In section 3, 
using some topological energy bounds on a compact domain, an asymptotic equation of state is 
derived. Then, in section 4, we make a conjecture about subleading terms in the EoS by exploiting 
some properties of generalised BPS Skyrme models. In section 5, the medium and large pressure 
regimes are checked numerically for the charge one, hedgehog configuration. Section 6 is devoted 
to a comparison with the equation of state (EoS) of the Walecka model. Einally, we summarise 
our findings. We shall frequently use the standard Skyrme field parametrization 

[/= cos^-|-isin^ n • a , n = (sinycos <I>, siny sin < 1 >, cos y) ( 1 . 6 ) 
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where X; ‘h are real field variables and a are the Pauli matrices. Further, we shall always assume 
that ?7 = 1, i.e., ^ = 0, is the vacuum value of the Skyrme field which it must approach at spatial 
infinity. 


II. AVERAGE PRESSURE AND AVERAGE CHEMICAL POTENTIAL 


The energy-momentum tensor in the BPS Skyrme model takes the form of a perfect fluid energy- 
momentum tensor and, therefore, it directly defines the field theoretical (microscopic) pressure 
and baryon chemical potential (i.e., the pressure and chemical potential densities). For a generic 
Skyrme Lagrangian, there is no such fluid description and one usually has a non-trivial energy- 
momentum tensor with different stresses in different directions. However, a pressure density 
p = \ T** and the corresponding average pressure P, i.e., the volume average of gT“ (summation 

assumed), as well as an average chemical potential fi may still be defined as proper thermodynamic 
(macroscopic) variables, i.e., obeying the required thermodynamical relations 




B 


-P, 


(II.l) 


dE 


= fj- 


V 


(II.2) 


Here E = E(V, P, B, fi) is the static energy of the Skyrme model, depending on the volume V, 
pressure P, baryon charge B and chemical potential fi. For this purpose, we first define the 
following generalized step function 


0(U) 


1 for U 
0 for U = 1 


(H.3) 


(for the standard parametrization (L6), Q{U) may be replaced by the standard step function 0(^)). 
Further, we define the locus function of a static skyrmion configuration U = Uq{x) as the pullback 
C/g(0(C/)) of 0(C/) under C/q, and the locus set of the skyrmion Uo{x), 


H = {x G I C/o*(0(U)) = 0([/o(x)) = 1}, 


(11.4) 


i.e., the set 0 C where the skyrmion is located (deviates from the vacuum). Now, let us consider 
the general Skyrme static energy functional 


E{V, P, B,fj,)= / d^x e[U] +P{ d^xQ{Uix)) -V] -fi{ / d^xBo - B 


(H.5) 


where e[U] is the energy density and P and are Lagrange multipliers. In particular, P imposes 
the condition that all possible solutions of the variational problem (H.5) must have volume V, i.e., 
f d^x 0(U(x)) = f^d^x = V. 0bviously, P and /i obey the thermodynamical relations (II.l) and 
(|H.2) by construction. We still have to show that the Lagrange multiplier P is indeed the average 
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pressure as defined above. For that we consider scaling transformations x* —>■ = (1 + A)x®. 


Then 5xU = Xx^diU. For the energy functional (II.5) we get off-shell, in hrst order in A, 

6xE = X [ Tiid^x - 3XP [ d^xQ{U{x)), 


(II.6) 


where we use the fact that the variation of the energy f d?x e[C/] under the scaling transformation 
gives the integral of the trace of the spatial part of the energy-momentum tensor m- But any 
solution of the Euler-Lagrange equations is a stationary point and, therefore, 5\E = 0 on-shell and 
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P = 




(11.7) 


which is exactly the average pressure definition. We remark that the lagrange multiplier imposing 
a fixed value for the volume does not uniquely determine the compact domain Q. In fact, using the 
volume preserving diffeomorphisms we get inhnitely many different sets H which result in different 
values of the trace integral (i.e., the pressure P). In other words, one has to fix the compact 
domain fl and, therefore, the skyrmion solution up to target space symmetries, to get a unique 
average pressure. In general, this is a very difficult task. However, for a spherically symmetric 
unit charge skyrmion, the obvious choice is such that the symmetry of the equilibrium solution is 
preserved while squeezed. Then, H is a three dimensional ball. Moreover, in the BPS limit the 
static energy functional is invariant under the volume preserving diffeomorphisms of the base space 
and therefore all SDiff related H give exactly the same pressure - as is the case for a perfect fluid. 


Let us also remark that the proper (thermodynamical) volume appearing in relation (II.l) is 


the geometrical volume of a topological soliton. It means that in a typical case such a volume 
tends to infinity in the equilibrium (zero pressure) limit since typical skyrmions are infinitely 
extended solutions. This leads to some difficulties if one wants to consider the mean-field energy 
density e = E/V at the equilibrium. Obviously, one gets e = 0 which is the same as for the vacuum 
configuration. Note, however, that compact skyrmions with finite volumes are known (for example, 
they are quite common in the BPS Skyrme model). 


III. ASYMPTOTIC EQUATION OF STATE 

Now we will use these definitions to compute thermodynamical functions in the asymptotic 
regime, i.e., in the limit of high pressure (density) or equivalently in the limit of small volume. In 
order to do that we need an expression that shows how the energy of skyrmions changes with volume 
and baryon density. This is a quite difficult problem as, in principle, it requires the solution of the 
full three dimensional problem. However, in the asymptotic regime we can use some previously 
found topological bounds. 

The relevant bounds for the quartic (Skyrme) and sextic terms are 


E, = ^ld^xTT > 3(27r2)4/3|^^ 


(HI.I) 




Eq= d^x TiLiLjLkY > — 


{ 111 . 2 ) 
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There is no bound for the quadratic part 

E2= f d^x Tv LiLi, 


(III.3) 


or for a general potential. Therefore, the asymptotic static energy for any Skyrme model can be 
approximated as (Ae = A^/(24)^ and Aq = 

d4/3 d2 d2 

E = X 2 E 2 + X 4 E 4 + X\^Ee + XqEo > A43(27r2)4/3_ + vr^A^— > vr^A^ —. (III.4) 

In the small volume limit, the quartic term provides a subleading contribution. The same holds 
for the potential term and the sigma model part, which affect the total energy even less. However, 
since in the small volume limit the sextic term governs the masses of skyrmions, one recovers the 
BPS limit. Hence, the bound is, in fact, saturated for sufficiently small volumes (high pressures). 
This is not the case for the quartic term. As we will see, the corresponding bound is not saturated. 
All this allows us to predict the asymptotic formula for the energy as 

,52 54/3 

with 


E = tt^A^-— + a 


+ V^O, 


(HI.5) 


a > 3(27r2)^/3A4. 


(HI.6) 


Then, the average pressure is 


D_ 4\2-2 , «-4/3 , /-4/3n 

P — TV X Pb + —Rb + o{Pb )) 


and the average baryon chemical potential 


O 4\2- , 4a_i/3 , /-1/3 n 

II = 2tt X PB + —Pb + oW ), 


(HI.7) 


(HI.8) 


where is the average baryon (particle) density. The average energy density is 

e = TT^X^p\ + ap^^'^ + o{p^q'^). (hi.9) 

In the leading approximation, we re-obtain the BPS Skyrme model asymptotic behaviour 


The resulting equation of state reads 


P = tt^X^pI, 

(HI.IO) 

A = ^TT^X^PB, 

(HI.II) 

£ = TT^X^p%. 

(HI.12) 

£ = P. 

(HI.13) 


Let us summarise the main findings; 
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1. The sextic term gives the main contribution in the high pressure limit. This means that this 
term should not be omitted if dense nuclear matter is considered. The asymptotic equation 
of state always has a universal (potential independent) form e = P. 

2 . The quartic, usual Skyrme term, gives a subleading contribution which modifies the equation 
of state at moderate pressures. The functional dependence is known, which is not the case 
for the multiplicative constant a, for which we have derived a lower bound. 

3. The potential and the sigma model part give contributions which are even subleading in 
comparison to the £^4 contribution. On the other hand, they may be significant close to 
nuclear saturation density. 

For a deeper insight into the role played by the sigma model part as well as by the potential 
term, we need to perform numerical computations. Nonetheless, one can get some help from the 
BPS Skyrme model which, as we will see below, can give some understanding on the functional 
mass-volume relation originating from the £2 and Eq terms. 


IV. THE BPS FLUID TOY MODELS 


Each part of the generalised Skyrme model has a specific number of spatial derivatives. This 
is a main - but clearly not the only - difference between them. As we have previously commented, 
the sextic term is special as it leads to a perfect fluid model. Here, we want to learn how a specific 
number of derivatives can change the equation of state, assuming that all terms are based on the 
same thermodynamical quantity - here the baryon density. Within this approach, we replace the 
sigma model term and the quartic Skyrme term by energy functionals based only on the baryon 
density, which is a natural quantity for the sextic term (Note, that similar non-linear models have 
been recently considered in [25] in the context of the baby Skyrme model in 2 -|- 1 dimensions). 

We start with the Skyrme (quartic) term 

£4 = A 4 TV {[L^, -^4 = (IV.l) 


which now is represented by a ’’four derivative” term (in the sense that under Derrick scaling 
— 7 - Ax^ it scales like A“^) constructed from the baryon density (current). Such a new term 
leads to the following energy density (for static solutions) 

e = (IV.2) 


and pressure 

P=1>^4pT, (IV.3) 

where we use that, for a general energy density based on baryon (particle) density, the pressure 
reads 


P = Pb 


ds 

dpB 


— e. 


(IV.4) 
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The corresponding equation of state reads 

P = (IV.5) 

which exactly coincides with the mean-field equation of state for the usual Skyrme term. Hence, 
the substitution completely reproduces the thermodynamical properties of the E 4 model. This 
leads to the conjecture that, at least to some extent, thermodynamical properties related to each 
term of the generalised Skyrme model can be similar to properties generated by the BPS fluid 
counterpart, i..e, by terms built out of the baryon density elevated to a certain power. Now, let us 
apply the same strategy to extract information about possible equations of state that may emerge 
from E 2 and Eq. In the case of the sigma model term we get 

Ts =-AaTr (L^L^) ^ 4 =(IV. 6 ) 

Thus, we find 

e = X2pT, (IV.7) 

P = -ypf > (IV.8) 

and the equation of state 

P = -ie. (IV.9) 

As expected, this part is responsible for a negative pressure, i.e., an attractive force (which is 
balanced by the inclusion of other stabilising terms). 

The issue of the potential is more subtle as it does not depend on the baryon density but, in 
general, on the Skyrme field. One finds that 

e = U, P = -U, (IV. 10) 

and 


P = -e. (IV.ll) 

However, for the step-function potential U = 0(t/), it only contributes as a numerical constant to 
thermodynamical quantities. Then, for example, the energy density is shifted by such a constant. 
Since in a first approximation, we can always model a potential by the step-function potential, we 
may expect that a constant density term shows up. 

Hence, finally we conjecture a theoretically motivated energy density-baryon density relation 
which can be valid not only for asymptotically high pressure but also in a medium pressure regime 

e = + j3 + j5 p^^^ + 0{l). (IV.12) 


This leads, for the charge one sector, to the following energy-volume formula 


1 


E = -|- a 


1 




+ /3V + /3 + 0(V). 


(IV.13) 


In the subsequent sections we will compare formula (IV.13) with numerical computations. One 
should keep in mind that, in contrast to the first two leading terms, which have been proven to 
emerge in the high pressure regime, the last two subleading terms are conjectured (or at best 
heuristically motivated) by analysing similar BPS models. 
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V. NUMERICAL RESULTS IN THE CHARGE ONE SECTOR 


The equation of state as well as the equivalence of the two pressure definitions can be tested by 
direct numerical computations in the charge one sector where the hedgehog ansatz can be assumed 

U = cos^ + i sin^n • (T, (V.l) 

where ^ = ^(r), n = (sin 0 cos </>, sin 0 sin (/), cos 0) is a unit vector covering the target subspace 
once, and a are the Pauli matrices. Then, the energy contributions are given by 

Eq = Att / r‘^dr2{l — cos^), (V.2) 


E 2 = 47r 


r‘^dr2 



(V.3) 


E4 = dvr 


r^dr 




(V.4) 


E(i=47r J r'^dr^sin^^Cr- (V.5) 

The aim is not only to verify the asymptotic formula for the EoS but also to get some insight 
into the role that is played by the quadratic and potential parts, for which we do not have any 
analytical (asymptotic) expression. For B = 1 skyrmions we solve a one-dimensional ODE with 
boundary conditions which guarantee the nontrivial topology and pressure 

^(r = 0) = vr, C{r = R) = 0. (V.6) 

In other words, we enclose the skyrmion in a hnite volume V = (47r/3)ii^. 

We use the collocation method [26j to determine the profile function ^ which minimises the 
Skyrme energy functional E with boundary conditions .^(0) = tt and ^{R) = 0. The associated 
Euler-Lagrange equations are solved on the interval Q < r < R. For the high pressure regime, we 
consider the interval 0.1 < ii < 1. We start with R = 0.1 as right boundary point and then increase 
in each collocation run the right boundary value by dr = 0.0001. The increment dr is chosen to 
be sufficiently small which allows us to check numerically the equivalence of field-theoretical and 
thermodynamical pressure. For the low pressure regime, we run simulations up to a maximal 
boundary point value of 10 with dr = 0.001. 

We use the nonlinear curve-htting function optimize.curve_fit from Python’s SciPy package m 
to fit the asymptotic energy and pressure formulas to our numerical data in the high pressure 
regime. 

The energy and length units have to be fixed by comparison with experimental nuclear physics 
data. In this article, we calibrate the usual Skyrme model without the sextic part using the 
approach of Refs. [2]. The energy and length units are tuned to match the experimental nucleon 
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(Mtv = 939 MeV) and delta (Ma = 1231 MeV) masses. This fixes the energy and length scale as 
follows 


127r^ X 


\ 1 

= 121 ? X 5.58MeV, 

2 

.4esky. 

_ _ 



= 0.755 fm, 


(V.7) 


and gives the well-known parameter values 

esky = 4.84, TV = 108 MeV, h = 46.8, and m-,^ = 138 MeV. (V.8) 


In this article, if not explicitly stated otherwise, the rescaled pion mass is set to 1. Note that 
the limitations of the parameter set (V.8) have been discussed previously |28H30j and different 
parameter choices have been explored in the literature, e.g. Refs. [SIITIEHIEI]. 

Our numerical results in section VB on BPS and near BPS Skyrme models are expressed in 


terms of the same energy and length units (V.7). However, note that for the full near-BPS Skyrme 


model T'o 246, discussed in section VB3 the coupling parameters of the BPS part are fixed by 
matching the BPS Skyrmion mass to one-fourth of the helium nucleus mass and the Skyrmion 
radius to the nucleon radius (see section VB3 and Ref. |32] for a detailed discussion). 


A. The perturbative Skyrme model 


We start with the usual, old Skyrme model without the sextic part. The reason for this is to 
find a clear signal from the quartic term which now gives the leading high pressure behavior and 
then understand the influence of the potential and sigma model terms as first subleading effects. 
The energy functional we minimise reads 

Eo 24 = A2TI2 -|- A4TI4 -|- XqEq, (V- 9 ) 


where the constants are 


A2 — 


247r2 ’ 


A4 


127r2’ 


Ao 


1 


(V.IO) 


or zero in the case where the pertinent term is omitted. Here we consider the usual Skyrme 
potential which is relevant for the conventional (perturbative) Skyrme model 


Eo= Tr(l-C/). 


(V.ll) 


For all possible submodels we first checked the equivalence between the average field-theoretical 
pressure and the thermodynamical pressure. We always found a perfect agreement between the 


pressures computed by Eqs. (H.l) and (H.7). It is also a precision test for our numerical compu¬ 
tations. 


1. Ei model 

We begin with the simplest submodel which contains only the quartic part. Such a limit of 
the Skyrme model has been considered before, however, from the instanton point of view m- 
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FIG. 1: (Color online) (a) Energy and (b) pressure of the B=1 skyrmion as a function of the volume for 
model. 


Obviously, due to the Derrick theorem there are no stable solitonic solutions at equilibrium (zero 
pressure), i.e., in the full infinite However, if we close the system in a hnite volume V, then 
skyrmions do exist. For a charge one soliton we have the following energy integral 


E4 


A 4 F /4 


127r2 


47r 




(V.12) 


This leads to a relatively simple equation of motion 


£ 


1 


cos ^ H—y sin ^ cos ^ = 0 


(V.13) 


dr"^ “ 

accompanied by the pertinent boundary conditions. Unfortunately, we were not able to solve it 
analytically. It is straightforward to notice that a solution on the segment [0, R] can be related to 
a solution on a different segment [0, AR] by the scale transformation r —)• Ar. Then, the energy 
scales with the factor A“^, which leads to 


= 1] ^ (V.14) 

where the coefficient a. can be understood as the energy in the unit volume. Hence, the pressure 
is a simple function of the volume 




3 U4/3 • 


(V.I5) 


Such an equation of state is confirmed by numerical computations. Indeed, we find that 


a = 1.853 or a = 924.4 MeV fm, 


(V.I6) 


where the first value is in Skyrme units while the second is in physical units. It is worth to compare 
this with the previously derived bound. Namely, 3(27r^)^/^/(I27r^) = 1.351. This shows that the 
bound is not saturated. The true solution is significantly above the bound. 
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Finally, we get an exact density-pressure equation of state which is valid for any value of the 
pressure 


e = 3P. 


(V.17) 


This result can be easily generalised to any topological sector and to any given compact fl. Indeed, 
let us consider the energy functional 

E4[f^] = [U,Ljf. (V.18) 

Obviously, it scales homogeneously under scaling transformations. Then again 


E4[E;fl] 


oiB,n 

1 / 1 / 3 ’ 


(V.19) 


where the constant a depends on the topological charge and a particular form of the compact 
manifold (assumed now to have unit volume). From bound pll.l ) we know that 


> 3(27r2)4/3B4/3. 


(V.20) 


However, as the energy functional is not invariant under the volume preserving transformations 
this constant should depend on the choice of H. In any case, it proves that the equation of state 
(V.17) is an off-shell universal equation of state for the pure quartic model. 

Finally let us remark that the bound is saturated if the geometries of the base and target spaces 
coincide. Indeed, if the base space is a three-sphere with a radius R then we get 


= 


—dvr — 
127r2 TR 


di; [ 2 sin^ 


(V.21) 


where ?// is the third angular coordinate on the sphere. This variational problem has a solution 
^ = Ip. Hence, the charge one skyrmionic solution is just the identity map between the base and 
target spheres. One can easily verify that in this case the bound is equal to the value of the last 
integral. This agrees with a classical result of Manton |34] where he shows that the identity map 
is a stable solution for the E24 Skyrme model for a sufficiently small radius R of the base space 
sphere. The geometrical reason why the M3 and §3 cases are so distinct is the following. The bound 
can be saturated if and only if the three eigenvalues of the corresponding strain tensor are 1) all 
equal and 2) constant which means that the corresponding map must be an isometry [23]. This is 
a very restrictive condition which, for the target space §3 ^ may be fulfilled only if the base space is 
also a three-sphere. 


2. Eq 4 model 

The pure quartic Skyrme model can be stabilised by the addition of a potential term. Then, the 
Derrick argument is evaded and an equilibrium solution exists. One peculiarity of this submodel 
(with the usual Skyrme potential) is that one gets a compact skyrmion, i.e., a topological soliton 










12 



xlO° 


1 , 

1 , 

1 , 

^ 0 , 

gO 

0 , 

0 , 

0 


- a/3=308.2,/3 = -13.6 


(a) 


6 8 

V[fm^ 

(b) 


10 12 14 16 




(c) (d) 

FIG. 2: (Color online) Energy and pressure of the B=1 skyrmion as a function of the volume for Eq 4 and 
E 24 model. 


which achieves its vacuum value at a finite distance. The size of the compact charge one skyrmion 
at equilibrium is 

R = 2.07 or i? = 1.563 fm, 

where the first number is in Skyrme units while the second in physical units. This sets the maximum 
volume (i.e., the volume at zero pressure) to Vmax = 15.9 fm^. 

The general asymptotic energy formula 

E = a^+pV^/^ + pV + o{V), (V.22) 

gives the following values for the parameters 

ao4 = 924.6 MeV fm, /3o4 = 13.4 MeV fm“^, /3o4 = 0. (V.23) 

In Fig.j^we plot the numerically computed energy (mass) and pressure as a function of the volume 
together with the theoretical formulas. We find a perfect agreement. We can also nse the theoretical 
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(a) (b) 




(c) (d) 

FIG. 3 : (Color online) Energy and pressure of the B =1 skyrmion as a function of the volume for i?024 model. 
Fig. (a)-(c) with the theoretical fit. Fig (d) with the best fitted curve. 


formula (V.22) to get the mass of the equilibrium skyrmion. Then, we get M 04 = 580 MeV which 
is only 3.5% above the true mass. Undoubtedly, the theoretical formula works very well for this 
submodel. 


3. E 24 model 

The conventional way for stabilising the E 4 model is to add the sigma model part. This is the 
usual massless Skyrme model. Now the equilibrium solution in the charge one sector is an infinitely 
extended skyrmion. Again, we want to fit the theoretical energy formula. The first issue is that 
we get a slightly smaller value of the leading term. Namely, 024 = 915.9 MeV fm. This means 
that we are still a bit away from the exact asymptotic regime. Apparently, one has to squeeze 
the skyrmion to higher pressures to reach the proper value for a. It shows that the sigma model 
part provides a significantly stronger attractive force than the potential term, which influences the 
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FIG. 4: (Color online) Comparison plots for the B=1 skyrmion: (a)-(d) mass and pressure as a function of 
the average baryon density; (e) mass as a function of pressure; (f)-(g) equation of state. 
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squeezed configurations and has a stronger impact on the energy-volume relation. 

In the fitting procedure we take this into account by imposing that a 24 = 924.1 MeV fm. Then, 

0:24 = 924.1 MeV fm, j324 = 0, ,024 = 311.5 MeV fm“^. (V.24) 

In Fig. we compare the numerical results and the theoretical curve. Again we find nice agree¬ 
ment in the high and medium pressure regime. For sufficiently small pressure (large volume) the 
theoretical formula is not valid any longer as the term gives a diverging contribution. 

4- Eq 24 model 

Finally, one can include all three terms. Then, after fitting the theoretical curve we get (with 
the same remark on the a constant which again does not take precisely its asymptotic value but 
a = 915.9 MeV) 

«024 = 924.1 MeV fm, /3o24 = —106.1 MeV fm“^, /3o24 = 458.0 MeV fm“^. (V.25) 

Note that the constant /?, which is related to the potential part of the action, takes a negative 
value. In fact, the potential and sigma model term contribute together to the subleading part of 
the energy-volume relation, which leads to this effect. It also suggests that the theoretical formula 
(the last two terms motivated by the BPS model analysis) is not quite adequate for this model. 
In fact, if we subtract the (corrected) leading term from the numerical energy then the best fit is 
found for V“ where a ~ 1/5 - see Fig. This leads to the following effective mass-volume formula 

^024 = + 7^1/5 + o(V^/®), (V.26) 

where 7 = 346.1MeV fm“^/®. This shows that in the full perturbative Skyrme model there is a 
rather strong mixing between the sigma model and potential part (and probable also with the quar- 
tic term) which effectively provides new subleading terms in the energy-volume relation. Therefore 
the theoretical curve proposed previously does not seem to work quite well in the full perturbative 
Skyrme model. 

In Fig. I^we plot numerical results describing how the energy (mass) and pressure depend on the 
average baryon density ns = 1/V. Furthermore, we plot the resulting equation of state relating 
the energy and pressure. 

It is also instructive to compare the mean-field energy density e with the true energy density 
distribution. This comparison makes no sense at the equilibrium (P = 0) where the soliton is 
infinitely extended. This means that the geometric volume is infinity V = 00 and the corresponding 
mean-field energy density is simply zero. However, for squeezed configuration the volume is always 
finite. Then, at least for the medium and high pressure regime the mean-field values obtained here 
for the charge one sector should coincide with their counterparts derived for a skyrmion on a torus 
(crystal). In Fig. [^we plot the local energy density for different values of the pressure together 
with the corresponding mean-field energy density. It is visible that the energy density has a global 
maximum whose value is significantly bigger than the average mean-field energy density. For small 



16 


pressure, this is not surprising, as the corresponding solution has a rather big volume. However, 
even for very large pressure {P = 1.432 x 10® MeV fm~®) the maximum of the energy density 
is approximately five times bigger that the average value. For medium pressure this difference 
increases. For example for P = 9.231 MeV fm“^ the maximum energy density is 20 times bigger 
than the average value. 

Let us finally remark that the leading term coefficient found in our numerical computations, 
a = 1.853 in Skyrme units (or a = 924.4 MeV in physical units), provides a bigger number (a 
better bound) not only than the analytical bound but also than the value found by Kutschera et. 
al. for a multi-skyrmion configuration being a cubic lattice of skyrmions, which in the Skyrme 
units is akut = 1-837 [35]. As in the asymptotic regime the multi-skyrmion should behave identical 
to the charge one skyrmion, we think that this constant gives the proper bound for high pressure. 


B. The near BPS Skyrme model 


In this subsection, we switch on the sextic term which, as we already know, governs the high 
density regime. The total energy has the following form 

Eo 246 = e (A 2 E 2 + X 4 E 4 + AqEo) -|- XqEq + XqEq. (V.27) 


The constants are 


A 2 — 


247r2 ’ 


A4 = 


1 

1^’ 


Aq — 


2„4 


127r2^ 


As = A^tt 


127r2’ 


Aq — 


127r2’ 


(V.28) 


where A and will be fixed later (or set to zero if the corresponding term is omitted). For the 
usual (perturbative) Skyrme model part, this corresponds to the same choice as before, but which 
the overall multiplication by the parameter e, which is assumed to have three values e = 1,0.1, 0.01. 
Again, the potential entering the perturbative Skyrme part is the standard potential (providing 
the physical mass for pions) 


Eo= TV(l-t/), 


while the BPS part of the potential is Eq = E^. 


(V.29) 


1. The BPS model - Eqq model 

In order to test the numerics we begin with the BPS Skyrme model, which can be solved 
analytically. The field equations of motion can be integrated to the following first order equation 

m 

X 127r^ 

^sin^^^r = -p^Eo + P-^, (V.30) 

where we assumed the hedgehog ansatz for the charge one soliton. Of course, the BPS equation 
can be found as the zero pressure condition [36|. The potential reads Eq = 4(1 — cos^)^. Using this 
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FIG. 5: (Color online) Local energy density distribution for the B=1 skyrmion in the E 024 model for different 
values of the pressure together with the corresponding mean-held energy density. 
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FIG. 6: (Color online) Energy and pressure of the B=1 skyrmions as a function of the volume for 
model. 
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equation, it is possible to compute the total energy integral and find the volume of the solution 
with a given pressure. Indeed, we find [H] 


E{P) 



sin^ ^ 


8{1 - cosC)^ + 

_ M 

^4(1-cos02 + P^’ 


(V.31) 


V{P) = 2Tr- [ dgsin^^^= ^ - . (V.32) 

^4(1 _ cos+ 

This gives the following mean-field equation of state 

, ^ p * ,2 /5 jOKMiiMidkiSi 

V2sO[{j,|,|}.{i,2}.-¥l/’ 

where pFq[{ai ,..., Op}, {6i,..., bg}, z] is a generalised hypergeometric function, and = l/37r^. 
Note, however, that the true energy density is not spatially constant and, in fact, differs quite 
significantly from its mean-field approximation. This is of some importance, e.g., for the application 
to neutron stars. For example, the mass-radius relation changes significantly if we switch from the 
true non-mean field energy density to the mean-field average energy density HU- 
The mean-field equation of state results, in the asymptotic region, in a simple energy-volume 
formula 



where 


E = a- + PooV + o{V), 


(V.34) 


a = 


1 

127r2 


A vr 


^ = 935.6 Mev fm^. 


/ 3 oo 


= 64.8 MeV fm"^. 

67r2 


(V.35) 


in Skyrme units and physical units (with the ANW calibration assumed). Here we chose A = 2 
and fi = 1. This asymptotical off-set can be compared with the equilibrium energy density 


£0 = ^ = 26.0 MeV fm-^. (V.36) 

37r4 

The non-zero value is due to the compact nature of the equilibrium solution. In fact, it has a finite 
size and a volume 


R = 



V = 2t? = 8.49 fm^ 


(V.37) 


Let us comment that for the most interesting full near-BPS Skyrme model we apply a different 
choice for parameters and different calibration scheme, which was previously used for neutron stars. 
In numerical computation we reproduce the mean-field equation of state for the BPS Skyrme model. 
For instance, we find 


a = 938.2 MeV fm^ and /3 = 61.5 Mev fm 


(V.38) 


which is in agreement with the analytical results within 0.3% and 5% respectively. 
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2. i ?026 model 

Now we consider a submodel for which the attractive force is provided entirely by the sextic 
term i.e., the quartic term is absent. Additionally we assume A = 2 and /r = 1. The leading term in 
the energy-volume relation is found in the numerical computations. Specifically, we reach a regime 
where 


de=o.oi — 938.4, d;£=o,i — 939.0, — 940.6. (V.39) 

Obviously, the leading part more correctly describes also a medium volume regime if we are closer 
to the BPS limit i.e., if we reduce the e parameter (see Fig. [^. We also get the subleading (linear 
in the volume) term. For the nearest BPS Skyrme considered here (e = 0.01) we find 

/3e=o.oi = 62.2 MeV fm-^ (V.40) 

with, however, a non-zero value for the off-set. It is also possible to fit -|- /dV curve. For 

e = 0.1 and e = 1 it leads to a negative value for /3, which seems to indicate that such a curve is 
rather not the right one. 


<?. Ao 246 model 


Finally, we consider the full near-BPS Skyrme model. As we are close to the BPS limit (for e = 
0.1 and e = 0.01) it is reasonable to apply a different choice for the calibration. Namely, we fit the 
parameters of the BPS part of the model to the mass of the helium nucleon mne/^ = 931.75 MeV 
and size of the nucleon tat = 1.25 fm. For the BPS Skyrme model we have 

Ebps = ( V - 41 ) 




R = 


V 2 7 


1/3 


2/i 7 


(V.42) 


Then, = 2317 Mev fm^, = 3372 MeV fm“^ in physical units or A^ = 8.14 and = 2.19 in 
the Skyrme units. In our numerical simulations we find the following values of the leading term in 
the energy-mass relation 


dte=i = 1929 MeV fm^, ae=o.i = 1908 MeV fm^, de=o,oi = 1905 MeV fm^, (V.43) 

which can be compared with the theoretical value a = = 1905 MeV fm^. 

As expected, the sextic term gives the leading behavior for the mass-volume formula in the 
asymptotic regime. The subleading contribution, emerging from the quartic part of the action, has 
the form with the following values for the constant Fig. 


«£=! = 924.7 MeV fm, Q;e=o.i = 92.7 MeV fm, ae=o.oi = 9.8 MeV fm, 


(V.44) 
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FIG. 7 : (Color online) Energy and pressure of B =1 skyrmion as a function of the volume for E0246 model. 
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FIG. 8: (Color online) Comparison plots for the B=1 skyrmion: (a)-(d) mass and pressure as a function of 
the average baryon density; (e)-(f) mass as a function of pressure; (g)-(h) equation of state. 
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which to a reasonable precision satisfies an obvious relation ag = e ae=i- This can be treated as a 
test for our numerics. 

In Fig. we compare the near-BPS Skyrme model (for three values of e) with the usual (per¬ 
turbative) Skyrme model. Of course, for a sufficiently large value of the average particle density 
the total energy (mass) as well as the pressure obey the relation which follows from the pure BPS 
model (sextic part), see Fig. [^(a), (c). For smaller values of ub (less squeezed configurations) i.e., 
closer to the equilibrium the relations are modified by a subleading behavior related to other terms 
of the action. In general, taking smaller e results in a flatter M — ns curve for small ns- This 
is perhaps related to the fact that the quartic term is more important for small densities and its 
strength is increased or reduced by a particular value of e. Furthermore, the M — hb and p — ub 
curves grow much faster for the full model (M ~ ub and p ~ n^) than in the perturbative Skyrme 
model (M ~ and p ~ n^B^). 


VI. COMPARISON TO THE WALECKA MODEL 


Let us now understand the results from the point of view of the Walecka effective model 
The pressure and energy density are given as 


.2 1 „2 1 


P= ^ 

2m2 *^47r2 


^4 


-tp - {m*ykF ] E*p + {my In - 

o / 


(VI.l) 


-L I 

2m2/^^2m2 " 47r2 


{2kl + {mykF) E*p - {my In ^ 


m^ 


where the in-medium Fermi energy (chemical potential) and in-medium nucleon mass are 


Ep = \l k% + (m*)2, rn* = m — -^n^. 

V ^ 


The baryon and scalar densities at T = 0 are 

m* 


PB = 


2k^ 


n. = 


TT^ 


kFEp — {m*y\n ^ ^ 


m^ 


(VI.2) 


(VI.3) 


(VI.4) 


It is straightforward to find the leading behaviour of the mean-field energy density as a function 
of the mean-field baryon density 


_ 1 5^ -2 , 3 /37r2y/^ _4/a X 


.4/3 ^ 


2my^ 4 V 2 ; ' 2 g2 ’ 


(VI.5) 


where the leading part comes from the u meson repulsion, the first subleading term is due to the 
free Fermion gas while the second subleading term comes from the scalar meson (more specifically, 
we use that in the asymptotic regime the scalar condensate tends to a constant value), 

n. = (VI.6) 


9o 


The results can be summarised as follows: 
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1. The asymptotic formulas for the equations of state in the full (i.e., containing the BPS part) 
Skyrme model and in the Walecka model coincide. Their functional dependence is exactly 
the same. The leading part implying e ~ P comes from the baryon density current. In the 
Skyrme model it enters via the sextic term while in the Walecka model it appears due to the 
oj meson repulsive interaction dominating at high density. So, we may detect the a; meson 
in the Skyrme action - as an emergent object hidden in the sextic term. 

2. Furthermore, there is a coincidence between the first subleading terms, which behave as 

In the Skyrme model it is generated by the usual Skyrme (quartic) part of the action while 
in the Walecka it is due to the ultra relativistic free fermion (nucleon) behaviour. 

3. Moreover, we got a constant term whose origin is found to be the potential (the Skyrme 
model) or scalar density condensation (the Walecka model). This gives further support to 
the idea that the potential corresponds to the cr meson and provides the attractive long 
distance force. 

4. Let us note that there are no further terms in the large density limit of the Walecka model 
which could correspond to the term generated by the sigma model part of the Skyrme model. 
Namely, a term behaving as with 4/3 > a > 0. 


The observed similarity between the equations of state can be used to get some quantitative 
insight into a; and a mesons within the Skyrme model framework, where such low energy particles 
are not explicitly included. If we compare the coefficients in the expression above we find 


TT 


,4x2 _ ^ ffg; 

2 m/, ’ 


(VI.7) 




(VI.8) 


I 


^ = 9 2 • 

2 gi 

The typical parameter values in the Walecka model are 


(VI.9) 


rriM = 939 MeV, = 783 MeV, ma = (500 — 600) MeV, 


(VI.IO) 


2 2 

^ = 14.717, ^ = 9.537. 

47r 47r 

This gives (after taking into account the Plank constant h = 197.3 MeV fm) 


(VI.ll) 


1 gl 

27r^ 


12 MeV fm^ 




362 MeV fm, 


1 

2 gl 


(120 - 175) MeV fm“3. 

(VI.12) 



25 


This can be compared with values derived for the generalized Skyrme model. For the BPS model 
(or near-BPS version with small e, for example equals to 0.01) }? and jj? for the pure BPS Skyrme 
model in the case of the step-function potential (with the parameter fit to the nuclear saturation 
density and the binding energy per nucleon of infinite nuclear matter) 

= 31 MeV W, /3 = 71 MeV fm“3. (VI.13) 

As we know, in the case of the usual (perturbative) Skyrme model, with the calibration assumed 
here one gets 

a = 924 MeV W. (VI. 14) 

Thus, for near-BPS Skyrme models this should be multiplied by the e constant, which leads to a 
rather small number in comparison to the Walecka model. 

VII. CONCLUSIONS 

In the present work we have analyzed thermodynamical properties of skyrmionic matter at zero 
temperature for the most general Skyrme theory that possesses both Poincare invariance and a 
standard hamiltonian. 

The main result is that the sextic term is responsible for the leading behavior at high pressure 
and therefore is unavoidable for any realistic application of the Skyrme model for a description of 
nuclear matter at high densities. As we found, the dependence of all global characteristics (total 
energy, pressure, average energy density) on the geometric volume tends to analytical relations 
derived for the BPS Skyrme model as the volume decreases. The same happens with the mean- 
field equation of state which approaches a much stiffer form, e = P, instead of e = 3P for the 
usual Skyrme model without the sextic part. Similarly, the mean-field energy density grows much 
more rapidly with the average baryon density - the previously found e ~ relation [35] is 
replaced by e ~ p^. This resolves a long standing discrepancy between the Skyrme model and 
other conventional models of nuclear matter. Indeed, the behavior generated by the sextic term 
coincides with the corresponding relation derived by Bethe and Johnson [39| (non-relativistic) as 
well as Walecka [37] (relativistic). Moreover, the inclusion of this term might not only modify the 
phase diagram of cold nuclear matter or neutron stars, but can also influence the scattering of 
nuclei m- It is also an intriguing observation that an effective low energy model of QCD, which 
here is represented by the Skyrme model, has an solvable and integrable limit (in the sense of the 
generalised integrability) at asymptotically high density (pressure). 

The first subleading contribution is provided by the Skyrme term, i.e., by the quartic part of 
the action, and it can modify the thermodynamical properties at medium densities (pressure). 
Although we analytically derived the functional dependence for this contribution, we found only 
a lower bound for the true value of the proportionality constant a. The true value was obtained 
numerically. The next subleading terms were found to be related to the attractive channels of the 
model, which emerge from the sigma model term and the potential. However, due to the strong 
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nonlinearity of the theory, there is a significant mixing in these subleading contributions which 
results in the appearance of some new, effective terms in the energy-volume relation. In particular, 
in the full £"0246 model, instead of terms and V - which are clearly visible in the £24 and £94 
submodels, we find an effective a G ( 0 , 1 ) dependence. 

Let us also note that, from the point of view of nuclear matter properties, our results are 
relevant for the high and medium pressure regimes. Dehnitely, the obtained mean-field equations 
of states are not applicable at the saturation density, where we got zero energy density. In fact, 
it is a common feature of Skyrme type models (with the exception of the BPS Skyrme model and 
the Eo 4 submodel) that £ = 1 solitons are infinitely extended at equilibrium (£ = 0 ). Then the 
equilibrium volume is inhnite, which leads to a vanishing mean-held energy density. This is not 
the expected behaviour for inhnite B. In the original Skyrme model, e.g, at inhnite B a crystal¬ 
like conhguration with a hnite volume per baryon number is conjectured to be the true minimizer 
even for small P. This does not mean that one cannot use hnite charge skyrmions to understand 
thermodynamical properties of nuclear matter or neutron stars. It simply means that one has to 
be very careful with the mean-held approximation. The full (non-mean held) densities are still hne. 
We want to underline that in our computations we solve the full held equations which provides 
us with full (spatially dependent) densities. Then they are averaged leading to the mean-held 
equation of state. This is an alternative approach to mean-held computations previously done for 
Skyrme-related models | 13 ] . 

If compared with another widely used effective model, i.e., the Walecka model, we found a 
correspondence between the cu meson interaction and the sextic term in the Skyrme model. Both 
generate the same mean-held equation of state, e = P. This again supports the crucial role played 
by the sextic term (and the BPS Skyrme model part) in the Skyrme framework. There is also 
a clear correspondence between the main subleading terms related to the ultra relativistic free 
fermion gas (Walecka) and the quartic Skyrme term (Skyrme model). The next subleading term in 
the mean-held equation of state in the Walecka model emerges due to the attractive scalar meson 
interaction. A similar constant term was found in the Skyrme model as well and it is generated 
by the potential. This provides some support to the natural expectation that the potential in 
the Skyrme model corresponds to the a meson in the Walecka model. However, the situation is 
much more involved as in the Skyrme model the term is accompanied by a term, which is 
generated by another source of the attractive force, namely the sigma model part. Furthermore, 
there is no such term in the Walecka model. In addition, the observed effective mixing between 
these attractive channels does not allow for a clear understanding of a possible emergence of the 
hidden scalar meson in the Skyrme model, in a similar fashion as it has been recently understood 
in the case of the emergent oj meson HD HU- In any case, it is an interesting feature of the Skyrme 
model that it does include some effects related to low lying mesons (here the u and a) although 
they are not introduced at the level of fields. 

There are several natural directions in which the present investigation can and should be further 
developed. First of all, one should find the mean-field equation of state in the limit of infinite nuclear 
matter. This can be achieved by considering the general Skyrme model, with the sextic term 
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included, on a torus. The main advantage of this set-up follows from the fact that the minimum 
energy solution occurs for a particular size (volume) of the torus. Therefore, the equilibrium 
solutions (at zero pressure) do not correspond to infinite geometrical volume which, as we already 
commented, is a drawback of finite charge solutions. As a consequence, equilibrium solutions 
possess finite mean-field energy and particle densities. Besides the derivation of the equation of 
state and understanding how it differs from the charge one results (which should be negligible in 
the high pressure limit), one can expect the appearance of some new phases (as, for example, a 
half-skyrmion phase). It would be desirable to fully discover the Skyrme phase diagram with the 
coupling constants as free parameters. 

Second, one can ask the more mathematical question what is the energy minimiser in a given 
topological sector on a given compact manifold with a certain volume (see for example [44 | ). 
Obviously, the form of the solution does not depend only on the volume of the manifold but 
also on its geometry (shape), although, as has been shown here, the asymptotic (small volume) 
behaviour is geometry independent. 

Third, it would also be desirable to extend our results to non-zero temperature, T > 0. The 
simplest way to accomplish this problem is to apply the caloron approach, i.e., to use euclidean 
periodic solutions in 4 dimension [35]. Alternatively, at least in the BPS limit, one can try to 
compute the partition function which is essentially given by an integration over the moduli space 
|46| . Unfortunately, for the BPS Skyrme model this group has infinite dimension which makes the 
problem quite involved. 


Acknowledgement 

The authors acknowledge financial support from the Ministry of Education, Culture, and Sports, 
Spain (Grant No. FPA2008-01177), the Xunta de Galicia (Grant No. INGITE09.296.035PR and 
Gonselleria de Educacion), the Spanish Consolider-Ingenio 2010 Programme CPAN (GSD2007- 
00042), and FEDER. AW thanks S. Krusch for discussion and D. Harland for inspiring remarks on 
the average pressure issue. Further, the authors thank J. Sanchez-Guillen for helpful comments. 


[1] T. H. R. Skyrme, Proc. Roy. Soc. bond. 260, 127 (1961); Nucl. Phys. 31, 556 (1962); J. Math. Phys. 
12, 1735 (1971). 

[2] G. S. Adkins, C. R. Nappi, and E. Witten, Nucl. Phys. B228, 552 (1983); G. S. Adkins and G. R. 
Nappi, Nucl. Phys. B233, 109 (1984) 

[3] G. J. Houghton, N. S. Manton, and P. M. Sutcliffe, Nucl. Phys. B510, 507 (1998); R. A. Battye and 
P. M. Sutcliffe, Nucl. Phys. B705, 384 (2005); Phys. Rev. C73, 055205 (2006); R. A. Battye, N. S. 
Manton, and P. M. Sutcliffe, Proc. Roy. Soc. bond. A463, 261 (2007); D. T. J. Feist, P. H. C. bau, 
and N. S. Manton, Phys. Rev. D87, 085034 (2013) 

[4] E. Braaten, b. Garson, Phys. Rev. bett. 56, 1897 (1986); Phys. Rev. D38, 3525 (1988) 

[5] b. Garson, Phys. Rev. bett. 66, 1406 (1991); Nucl. Phys. A535, 479 (1991); T. S. Walhout, Nucl. Phys. 
A531, 596 (1991) 



28 


[6] O. V. Manko, N. S. Manton, S. W. Wood, Phys. Rev. C76, 055203 (2007); R. A. Battye, N. S. Manton, 
P. M. Sutcliffe, S. W. Wood, Phys. Rev. C80, 034323 (2009) 

[7] P. H. C. Lau, N. S. Manton, Phys. Rev. Lett. 113 (2014) 232503 

[8] I. Klebanov, Nucl. Phys. B262, 133 (1985); L. Castillejo, P. S. J. Jones, A. D. Jackson, J. J. M. 
Verbaarschot, and A. Jackson, Nucl. Phys. A501, 801 (1989); M. Kugler and S. Shtrikman, Phys. 
Lett. B208, 491 (1988); M. Kugler and S. Shtrikman, Phys. Rev. D40, 3421 (1989); J. M. Speight, 
Commun. Math. Phys. 332 (2014) 1, 355. 

[9] C. Adam, J. Sanchez-Guillen, A. Wereszczynski, Phys. Lett. B691, 105 (2010); C. Adam, J. Sanchez- 
Guillen, A. Wereszczynski, Phys. Rev. D82, 085015 (2010). 

[10] C. Adam, G. Naya, J. Sanchez-Guillen, J. M. Speight, A. Wereszczynski, Phys. Rev. D90 (2014) 045003 

[11] C. Adam, C. Naya, J. Sanchez-Guillen, R. Vazquez, A. Wereszczynski, Phys. Rev. C92 (2015) 025802 

[12] C. Adam, C. Naya, J. Sanchez-Guillen, A. Wereszczynski, Phys. Rev. Lett. Ill (2013) 232501; Phys. 
Rev. C88 (2013) 054313. 

[13] A. Jackson, A. D. Jackson, A. S. Goldhaber, G.E. Brown, L. G. Castillejo, Phys. Lett. B154 (1985) 
101 

[14] U.-G. Meissner, 1. Zahed, Z. Phys. A327 (1987) 5 

[15] R. Alkofer, H. Reinhardt, H. Weigel, Phys. Rept. 265 (1996) 139 

[16] W.-G. Paeng, H. K. Lee, M. Rho, C. Sasaki, Phys. Rev. D88 (2013) 105019; M. Rho, Int. J. Mod. 
Phys. A25 (2010) 5040. 

[17] J.-H. Jung, U. Yakhshiev, H.-C. Kim, P. Schweitzer, Phys. Rev. D89, 114021 (2014) 

[18] G.-J. Ding, M.-L. Yan, Phys. Rev. C75 (2007) 034004 

[19] E. Bonenfant, L. Marleau, Phys. Rev. D82 (2010) 054023; E. Bonenfant, L. Harbour, L. Marleau, 
Phys.Rev. D85 (2012) 114045; M.-O. Beaudoin, L. Marleau, Nucl. Phys. B883 (2014) 328 

[20] J. M. Speight, J. Geom. Phys. 92 (2015) 30 

[21] M. Gillard, D. Harland, J. M. Speight, Nucl. Phys. B895 (2015) 272 

[22] S. B. Gudnason, M. Nitta, Phys. Rev. D91 (2015) 085040; S. B. Gudnason, M. Nitta, Phys. Rev. D91 
(2015) 045018; S. B. Gudnason, M. Nitta, Phys. Rev. D91 (2015) 045027; S. B. Gudnason, M. Nitta, 
Phys. Rev. D90 (2014) 085007; S. B. Gudnason, M. Nitta, Phys. Rev. D89 (2014) 025012 

[23] D. Harland, Phys. Lett. B728 (2014) 518 

[24] C. Adam, A. Wereszczynski, Phys. Rev. D89 (2014) 065010 

[25] J. Ashcroft, M. Haberichter, S. Krusch, Phys. Rev. D91 (2015) 105032 

[26] U. Ascher, J. Christiansen, R. Russell, ACM TOMS, 7 (1981) 209; U. Ascher, J. Christiansen, R. 
Russell, Math, of Comp., 33 (1979) 639 

[27] “Scipy, scientific tools for python, version 0.6.0”, http://www.scipy.org 

[28] R. A. Battye, S. Krusch, P. M. Sutcliffe, Phys. Lett. B626 (2005) 120 

[29] C. Houghton, S. Magee, Phys. Lett. B632 (2006) 593 

[30] R. A. Battye, M. Haberichter, S. Krusch, Phys. Rev. D90 (2014) 125035 

[31] D. Foster, N. S. Manton, Nucl. Phys. B899 (2015) 513. 

[32] C. Adam, C. Naya, J. Sanchez-Guillen, R. Vazquez, A. Wereszczynski, Phys. Lett. B742 (2015) 136 

[33] J. M. Speight, Phys. Lett. B659 (2008) 429 

[34] N. S. Manton, Commun. Math. Phys. Ill, 469 (1987) 

[35] M. Kutschera, C. J. Pethick, D. G. Ravenhall, Phys. Rev. Lett. 53 (1984) 1041 

[36] D. Bazeia, L. Losano, R. Menezes, J. C. R. E. Oliveira, Eur. Phys. J. C51, 953 (2007); D. Bazeia,L. 
Losano, R. Menezes, Phys. Lett. B668, 246 (2008); A. N. Atmaja, H. S. Ramadhan, E. da Hora, 
arXiv:1505.01241, 

[37] J. D. Walecka, Ann. Phys. 83 (1974) 491; B. D. Serot, J. D. Walecka, Int. J. Mod. Phys. E6 (1997) 



29 


515. 

[38] A. Schmitt, Dense Matter in Compact Stars, Lecture Notes in Physics 811, Springer Verlag, Heidel¬ 
berg 2010. 

[39] H. A. Bethe, M. Johnson, Nucl. Phys. A230 (1974) 1 

[40] D. Foster, S. Krusch, Nucl. Phys. B897 (2015) 697 

[41] C. Adam, T. Klahn, C. Naya, J. Sanchez-Guillen, R. Vazquez, A. Wereszczynski, Phys. Rev. D91 
(2015) 125037 

[42] D. Foster, P. M. Sutcliffe, Phys. Rev. D79, 125026 (2009). 

[43] G. Kalbermann, Nucl. Phys. A612 (197) 359; R. Ouyed, M. Butler, Astrophys. J. 522 (1999) 453; P. 
Jaikumar, R. Ouyed, Astrophys. J. 639 (2006) 354. 

[44] N. S. Manton, P. J. Ruback, Phys. Lett. B181 (1986) 137 

[45] M. Nowak and 1. Zahed, Phys. Lett. B230, 108 (1989); K.J. Eskola and K. Kajantie, Z. Phys. C44, 
347 (1989); Y. Shnir, G. Zhilin, Phys. Lett. B723 (2013) 236. 

[46] N. Manton, Nucl. Phys. B400 (1993) 624. 



